椭圆轨道卫星对地全球覆盖电推进控制

针对椭圆轨道卫星近/远地点的星下点对全球或特定纬度区域的访问问题，提出一种连续小推力下的对地覆盖控制策略。首先，推导了自然摄动对卫星拱线变化的影响，并探讨了进行小推力覆盖控制的必要性。然后，针对燃料消耗的优化问题，将控制方程展开成含傅里叶级数的形式，用以获得便于星上计算的解析形式的次优解，同时探讨了截取阶数与优化程度的关系。在进行拱线控制的同时，通过合理设置约束，对椭圆轨道的近地点高度进行保护，确保卫星安全运行。仿真结果表明，提出的方法能够以适当的燃料消耗代价实现椭圆轨道的近/远地点的全球覆盖控制或特定纬度区域的反复推扫，且控制力在可接受的范围内。     

同步双小行星系统共振轨道设计

研究了同步双小行星系统中共振轨道的设计方法及演化规律。首先，基于双椭球模型建立探测器运动方程，并给出共振轨道初值选取方法。然后，利用改进并行打靶法，提出一种双小行星系统平面共振轨道两步修正方法。同时结合稳定性理论及分岔理论，给出双小行星系统三维共振轨道生成和延拓方法；最后，以双小行星系统1999KW4为例，设计了共振比为1∶1,1∶2,1∶3,1∶4,2∶3的平面和空间共振轨道族，并分析了共振轨道的特性及轨道周期和轨道能量的变化规律。给出的双小行星系统中共振轨道的设计方法具有普适性，对未来双小行星系统探测任务中的轨道设计具有一定的参考意义与借鉴价值。     

Spectroscopy peculiarities of thermal plasma of electric arc discharge between electrodes with Zn admixtures

Plasma of the free burning electric arc between Ag–SnO₂–ZnO composite electrodes as well as brass electrodes were investigated. The plasma temperature distributions were obtained by Boltzmann plot method involving Cu I, Ag I or Zn I spectral line emissions. The electron density distributions were obtained from the width and from absolute intensity of spectral lines. The laser absorption spectroscopy was used for measurement of copper atom concentration in plasma. Plasma equilibrium composition was calculated using two independent groups of experimental values (temperature and copper atom concentration, temperature and electron density). It was found that plasma of the free burning electric arc between brass electrodes is in local thermodynamical equilibrium. The experimental verification of the spectroscopic data of Zn I spectral lines was carried out.     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

Families of halo orbits in the elliptic restricted three-body problem for a solar sail with reflectivity control devices

Solar sail halo orbits designed in the Sun-Earth circular restricted three-body problem (CR3BP) provide inefficient reference orbits for station-keeping since the disturbance due to the eccentricity of the Earth’s orbit has to be compensated for. This paper presents a strategy to compute families of halo orbits around the collinear artificial equilibrium points in the Sun-Earth elliptic restricted three-body problem (ER3BP) for a solar sail with reflectivity control devices (RCDs). In this non-autonomous model, periodic halo orbits only exist when their periods are equal to integer multiples of one year. Here multi-revolution halo orbits with periods equal to integer multiples of one year are constructed in the CR3BP and then used as seeds to numerically continue the halo orbits in the ER3BP. The linear stability of the orbits is analyzed which shows that the in-plane motion is unstable while the out-of-plane motion is neutrally stable and a bifurcation is identified. Finally, station-keeping is performed which shows that a reference orbit designed in the ER3BP is significantly more efficient than that designed in the CR3BP, while the addition of RCDs improve station-keeping performance and robustness to uncertainty in the sail lightness number.     

地球同步轨道薄膜太阳帆的姿态轨道控制方法

薄膜太阳帆（FSS）是集推进、发电和姿轨控功能于一体化的超大型挠性太阳帆式航天器，通过调整薄膜反射率产生可变推力和力矩，实现其姿态和轨道运动控制。结合薄膜太阳帆在地球同步轨道运行时的受力特性进行了轨道漂移分析。通过建立薄膜太阳帆动力学模型及受力模型，提出了调整帆面角度轨道修正方法以及基于薄膜光压力矩角动量卸载的长期在轨对日定向面内双轴动量轮稳定控制方法。通过系统仿真验证表明所提的轨道修正和对日定向控制方法是合理有效的，可使薄膜太阳帆长期在定点位置维持对日定向。     

基于TLE的弹道系数计算方法及应用分析

针对仅使用两行要素（Two Line Element，TLE）作为数据源的应用需求，研究了基于TLE轨道衰减的弹道系数计算方法。介绍了一种常用的基于两组TLE的直接计算法，分析TLE选取间隔对结果精度的影响；提出了一种基于多组TLE的迭代计算方法，以降低异常TLE对计算结果的影响；从弹道系数计算效果、在再入预报中的应用等方面对这两种方法进行比较分析。结果表明，两种方法各有优劣，基于多组TLE的迭代计算法稳定性更高、受TLE精度的影响更小；由于数据区间更短，基于两组TLE的计算结果对短期轨道衰减特性反应得更准确，用于临近再入时的预报效果更好。     

An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincare′sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L_1, L_2, and L_3 points as well as Lissajous orbits for L_1 and L_2 points in the EarthMoon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.     

一种低轨双星雷达信号无模糊频差估计算法

针对低轨双星组合接收低脉冲重复频率(LPRF)雷达信号的频差(DFO)估计存在模糊的问题，提出一种基于时差差分的频差无模糊高精度估计算法，该算法通过对高精度的时差测量值进行中心差分获得无模糊的频差估值。仿真结果表明，新算法能对PRF低至300 Hz的雷达信号进行无模糊的频差估计，估计精度在信噪比高于15 dB时逼近克拉美劳下界，可有效扩大低轨双星时/频差定位系统对LPRF信号的适用范围。     

Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves

We have used the radio occultation (RO) satellite data CHAMP/GPS (Challenging Minisatellite Payload/Global Positioning System) for studying the ionosphere of the Earth. A method for deriving the parameters of ionospheric structures is based upon an analysis of the RO signal variations in the phase path and intensity. This method allows one to estimate the spatial displacement of a plasma layer with respect to the ray perigee, and to determine the layer inclination and height correction values. In this paper, we focus on the case study of inclined sporadic E (E_s) layers in the high-latitude ionosphere based on available CHAMP RO data. Assuming that the internal gravity waves (IGWs) with the phase-fronts parallel to the ionization layer surfaces are responsible for the tilt angles of sporadic plasma layers, we have developed a new technique for determining the parameters of IGWs linked with the inclined E_s structures. A small-scale internal wave may be modulating initially horizontal E_s layer in height and causing a direction of the plasma density gradient to be rotated and aligned with that of the wave propagation vector k. The results of determination of the intrinsic wave frequency and period, vertical and horizontal wavelengths, intrinsic vertical and horizontal phase speeds, and other characteristics of IGWs under study are presented and discussed.